Towards photostatistics from photon-number discriminating detectors 
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We study the properties of a photodetector that has a number-resolving capability. In the absence 
of dark counts, due to its finite quantum efficiency, photodetection with such a detector can only 
eliminate the possibility that the incident field corresponds to a number of photons less than the 
detected photon number. We show that such a non-photon number-discriminating detector, however, 
provides a useful tool in the reconstruction of the photon number distribution of the incident field 
even in the presence of dark counts. 
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I. THE PROBLEM 

With the recent advent of linear optical quantum com- 
putation (LOQC), interest in photon number-resolving 
detectors has been growing. The ability to discrimi- 
nate the number of incoming photons plays an essen- 
tial role in the realization of nonlinear quantum gates in 
LOQC Jj Jp lp i B as well as in quantum state prepara- 
tion [E Eft & 13 ■ Such detectors are used to post-select 
particular quantum states of a superposition and con- 
sequently produce the desired nonlinear interactions of 
LOQC. In addition, the most probing attacks an eaves- 
dropper can launch against a typical quantum cryp- 
tography s yste m exploit photon-number resolving capa- 
bility [!(], [HI- Recent efforts in the development of 
such photon number-resolving detectors include the vis- 
ible hghtphoton counter [Ta Il3[ . fiber- loop detectors 
[bi. HH llfl lril Il8l ] , and superconducting transition edge 
sensors Standard photodetectors can measure only 
the presence or absence of light (single-photon sensitiv- 
ity), and generally do not have the capability of discrim- 
inating the number of incoming photons (single-photon 
number resolution). There have been suggestions of ac- 
complishing single-photon resolution using many singlc- 
photon-sensitive detectors arran ged in a detector array 
or detector cascade 0, For example, the 

VLPC (visible light photon counter) is based on a con- 
fined avalanche breakdown in a small portion of the total 
detection area and, hence, can be modeled as a detector 
cascade 0. Then again, fiber-loop detectors may be 
regarded as a detector cascade in the time domain. 

Let us suppose a number-resolving detector detects two 
photons in a given time interval. If the quantum effi- 
ciency of the detector is one, we can be certain that two 
photons came from the incident light during that time 
interval. If the quantum efficiency is, say, 0.2, what can 
we say about the incident light? In the absence of dark 
counts, we can safely say only that there were at least 
two photons in the incident beam. In this case, photode- 
tection rather has a non-photon number discriminating 



feature, since the conditional probability of not having 
zero, or one photon in the incident pulse is zero. Can we 
say more than that? For example, what is the probability 
that the incident pulse actually corresponds to two pho- 
tons, or three? In this paper we attempt to answer these 
questions and discuss the effect of quantum efficiency on 
the photon-number resolving capability. 

The semiclassical treatment of photon counting statis- 
tics was first derived by Mandel [25j . He showed that the 
photon counting distribution is Poissonian if the intensity 
of the incident light beam is constant in time. The full 
quantum mechanical description has been formulated by 
Kelley and Kleiner [2(|, and Glauber [2l}. Here, we start 
from a formula based on the photon-number state rep- 
resentation, which was first derived by Scully and Lamb 
|28l | . For a given quantum efficiency of the photodetector, 
the photon counting distribution can be written as 

m = jt( n X k (l-vT- k S(n), (1) 

n=k ^ ' 

where P{k) is the probability of detecting k photons, r\ is 
the quantum efficiency, and S(n) is the probability that 
the source (incident light) corresponded to n photons. It 
is assumed that the quantum efficiency r\ does not depend 
on the intensity of light. This distribution then leads to a 
natural definition of P(k\n) as the conditional probability 
of detecting k photons given that n photons being in the 
source as 

P(k\n)= (fyv k a-v) n - k - (2) 

Contrariwise, suppose we detect k photons. If the de- 
tector has a quantum efficiency much less than unity, very 
little can be said about the actual photon number associ- 
ated with the source. We may, however, deduce that the 
source actually corresponded to n photons with a cer- 
tain probability. This probability should be high if the 
quantum efficiency is high. In the following, we will use 
Bayes's theorem in order to extract information about 
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the source from the conditional photon-counting statis- 
tics. Armed with the conditional probability P(k\n), we 
may write 



Q(n\k) 



P(k\n)S(n) P(k\n)S(n) 



P(k) 



(3) 



Now we will identify the conditional probability, 
Q(n\k) with the probability that an ideal detector would 
have detected n photons, given that k photons are de- 
tected by the imperfect detector. It is given as the prod- 
uct of the probability of the source corresponding to n 
photons before the measurement and the probability of k- 
photon detection, given that n photons are in the source, 
normalized to the probability of having k photons de- 
tected. 

It is evident that generally we cannot find Q(n\k) if 
we do not know S'(n)-the distribution of the source. The 
exceptions are: (i) If n < k, then P(k\n) = 0. Therefore, 
Q(n\k) = for n < k no matter what S(n) is. (ii) If the 
quantum efficiency is unity, that is rj = 1, then P(k\n) = 
when n ^ k, and P(k\n = k) = 1. Since P(k) = S(n = 
k), it leads to Q(n = k\k) = 1 and Q(n ^ k\k) = 0. 
Hence, without knowing S(n), we only have information 
about Q(n\k) when n < k or 77 = 1. Such a detector, 
then, readily excludes certain numbers of photons in the 
source rather than determining them. Without a prior 
knowledge of the source, if n ^= 1, we call this a non- 
photon number- discriminating detector. 

Obviously, the difficulty in talking about this condi- 
tional probability lies in the fact that the photodetection 
is destructive. We will elaborate this subtlety using a 
hcraldcd-photon setup depicted in Fig. ^ The unitary 
operator Ur transforms the input state in the following 
way: 



a„|n)i|n) 2 . 



(4) 



n=0 



n=0 



Suppose now the photon number in the mode 1 is mea- 
sured by the perfect detector and the photon number in 
the mode 2 is measured by the imperfect detector. It 
is now simple to understand the conditional probability 

perfect detector 

V 




imperfect detector 



FIG. 1: Schematic of the heralded-photon setup. 



of finding n photons in the mode 1, given that k pho- 
tons are detected in the mode 2. In this way, we can 
have an operational definition of the conditional prob- 
ability Q(n\k). It is easy to show that the conditional 
probability of finding n photons in mode 1 given that k 
photons are detected in mode 2 is given by precisely the 
same expression for the quantity Q(n, k) as Eq. (3). [Note 
that the transformation given by Uh is far from quantum 
cloning, since the result of U is far from being equal to the 
cloning transformation (J2^=o a n\n)i) ® (X)^=o a n\n)2)- 



II. POVM DESCRIPTION 

Suppose an initial quantum state \ipi n ) of incident light 
is represented in the photon number states, as 



\lpi; 



J2a n \n). 



(5) 



n=0 



The probability of finding n photons is then given by 
S(n) = I (n\ip) | 2 = |a n | 2 . We may describe the detection 
process of the perfect detector by entangling the source 
with the detector ancilla |0) 



where 



\i>f) = U\ip in )\Q) =^2a n \n) Q \a n ) M , 



(6) 



(7) 



is the unitary operator implementing the von Neumann 
measurement, and we identified the quantum (Q) and 
measurement device (M) Hilbcrt spaces with appropriate 
subscripts. 

The probability to detect k photons, P(k), can now be 
calculated as 

P(k) - Tr Q (\a k ) M (a k \\^ f )(^ fl \) = \a k \ 2 , (8) 

as expected. Of course, the same result would have been 
obtained as P(k) = Tr(MkPi n ) using a Positive Operator 
Valued Operator Measure (POVM) 



M k = \k)(k\ 



(9) 



The present unitary treatment, however, is more trans- 
parent and allows a straightforward extension to imper- 
fect detectors. 

For an imperfect photodetector with a finite quantum 
efficiency, the POVM can be written as 



M k = J2 p ( k \n)K)( On , 



(10) 



n—k 



with P(k\n) given by Eq. We can verify this by 

calculating as before 



H' f ) = u\ij in )\o)M\o)E , 
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except that the entanglement operator U now acts on the 
joint Hilbert space of the quantum system, a measure- 
ment device, and an environment \0)e, which is where 
our lost photons go. Thus, the unitary operator can be 
given by 



U = 



where ip is an arbitrary phase and is the creation op- 
erator for the relevant mode of the environment. Similar 
to Eq. iJSJ, we now have 



P(k) = Tiqe (\a k ) M (ak\\ipfi)(^fi\) 



= £KI 2 



^(1-7?) 



n — k 



(12) 



supporting our Eq . UlUfl as the correct POVM. Note that 
any relative phase information between the photon num- 
ber states disappeared. 

The POVM effectively transfers the initial quantum 
state of light into a certain quantum state (as a result of 
fc-photon detection) by the transformation 



Pf-k — MkPin 



(13) 



where pi„ = | - 0m)(V'm|- The probability that k photons 
are detected can then also be written as 



P{k) = Tr( P/ _ fc ) 

= Tr(M fc |V> in ><^„|) 
= £i>(ft|n)S(n). 



(14) 



On the other hand, the joint probability that the source 
state is \n) and k photons are detected can be expressed 
by the probability 



P(n,k) = Tr E ((n,ak\tpf}(i>f\n,a k }) 
= (n|p/_ fc |n) = P{k\n)S(n) 

which defines the quantity Q(n\k) as 

{a n \p^k\a n ) 
1 ' ' Tr( P/ _ fc ) 



(15) 



(16) 



In particular, the quantity Q(n = k\k) corresponds 
to the confidence of the state preparation in a situation 
that we have more than two modes for the input and the 
photon number in one of the modes is measured |22j | . 



III. AN EXAMPLE 



A coherent state \a) is expressed in terms of number 



state as \a) 



72 ■ 



£ n (a"/V<)|n) |29j. The prob- 
ability of finding n photons in \a) is then given by a 
Poisson distribution: 



S(n) = \(n\a)f 



(17) 



where n — \a\ 2 . Now the probability of detection of k 
photons with an imperfect detector is written as |28l | 



P(k) = J2 



n—k 



i-k ( 



v k (i- v y 

k I TV. 



Efl + k\ n l n k 

1=0 



k\ 



k J(l + k)\ 



- fi (ny) k m . 



k\ 



{nri) h 
k\ 



-nr; 



(18) 



That is, the photon counting statistics are simply a Pois- 
son distribution where the average number of photons 
detected is the average number of incident photons mul- 
tiplied by the quantum efficiency of the detector-given 
by k = nr\ as expected. 

Now let us consider the conditional probability Q{n\k). 
Using Eq. ©, it is simply found as 



Q{n\k) 



(l)rj k (l - i 1 ) n - k e- n n n /n\ 
e- n ^{nrf) k /k\ 



{n-k)\ 



l\ 



-I 



(19) 



where / = n — k and / = ?i(l —77). Therefore, the form 
of Q(n\k) is again a Poisson distribution as a function of 




Let us take a simple example of an input coherent state 
in order to see how the value of Q(n\k) behaves. In doing 
this, we assume that the dark count rate is small enough 
that we can ignore it, and the quantum efficiency is in- 
dependent of the number of incident photons. 



FIG. 2: The conditional probability Q(n\k = 1) for a coherent 
state with n = 1 as a function of n. (a) For detector efficiency 
n = 0.2, Q(l|l) = 0.45, Q(2|l) = 0.36. (b) For detector 
efficiency 77 = 0.9, Q(l|l) = 0.90, Q(2|l) = 0.09. 
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FIG. 3: The conditional probability Q(n\k = 1) for a coher- 
ent state with n = 0.2 as a function of n. (a) For detector 
efficiency 77 = 0.2, Q(l|l) = 0.85, Q(2|l) = 0.14. (b) For 
detector efficiency 77 = 0.9, Q(l|l) = 0.98, Q(2\l) = 0.02. 



n — k with the average number of photons that fail to 
register given as n(l — 77). 

Figure shows the behavior of Q(n\k) as a function of 
n. for the case of coherent state input with n = 1, when 
one photon is detected. We can see that if the quantum 
efficiency of the detector is 0.2, with the detection of one 
photon, then the chances are more than 50% that the 
input contains more than two photons (Q(l|l) = 0.45). 
Obviously, Q(l|l) increases for higher values of the quan- 
tum efficiency (Q(l|l) = 0.9 for 7/ = 0.9). 

On the other hand, for a given quantum efficiency of 
the detector, Q(l|l) gets large as the average number of 
photons in the input decreases. For n = 0.2, Q(l|l) = 
0.85 even if the detector efficiency is 0.2, as depicted in 
Fig- El I 11 this case, with the detection of one photon the 
chances are 85% that the detection result corresponds to 
the correct value. 



IV. RECONSTRUCTION OF 
PHOTON-NUMBER DISTRIBUTION 

So far we have seen that we cannot say much about 
Q(n\k) if we do not know at least the form of S(n). In 
other words, using an imperfect detector, even after a 
detection of k photons, we do not gain much informa- 
tion about the number of incident photons without a 
prior knowledge of the incident field. Hence, the problem 
comes back to reconstruction of photon number distribu- 
tion of the incident field, S(n), from the photon counting 
probability, P(k) by a sequence of identical inputs and 
measurements. For this purpose, we may write Eq. (J2J 
in a matrix form as 

P = PS, (20) 



where 



P(l) 

P(2) 


, s = 


f 5(1) 
S(2) 


V ■ J 




{ ■ J 



and each component of the matrix P represents the con- 
ditional probability P(k\n) as 

P kn ^P(k\n)= (fj v k (l- V ) n - k . (22) 

Then the inverse of P can be found explicitly; it is given 
by 

(P-% m =(™)r 1 - m (v-ir- n , (23) 

which, in the absence of dark counts, determines the 
photon-number distribution of the source from the 
photon-counting distribution. Note, however, that for 
a fixed n there are alternating signs in (P _1 ) mn as the 
index m changes, so that this inversion formula can yield 
unphysical results depending on the probability distri- 
bution P{k). Put another way, for small efficiencies the 
matrix P is ill-conditioned and P 1 will contain matrix 
elements that can get very large in absolute value and 
can have either sign. This means that even if the cal- 
culated S(n)'s are all non- negative, a good inversion is 
still not guaranteed, since the inversion formula i|23|) is 
highly sensitive to small changes of P(k), amplifying the 
statistical and other noise which is inevitably present in 
any actual measurement. The fact that the inverse ma- 
trix P _1 is also a upper triangular matrix implies that 
S(n) is determined by P(k) with k > n. In other words, 
the tail of the measurement result, P(k), where the sta- 
tistical noise can be expected to be greatest, affects the 
inversion significantly. The validity conditions for the 
inversion formula shall be discussed elsewhere. 

Instead, we can carry out the following recipe: 1) Seek 
an underlying distribution S for which PS matches the 
observed frequencies according to a x 2 goodness-of-fit 
test. 2) Then, of all such distributions, seek the one 
with maximum entropy; i.e., pick the most likely under- 
lying distribution, which is statistically consistent with 
the data. Such a strategy can be implemented as a ge- 
netic algorithm, which we will discuss in a forthcoming 
paper [3fl |. Using this method, the reconstruction of S 
with a data set from an actual experiment is depicted 
in Fig. 0] The plot in Fig. [3Ji shows the probability dis- 
tribution of the detected photons. The raw count data 
was normalized by the total number of counts (data pro- 
vided courtesy of NIST 0. The result of reconstruction 
is shown by the plot Fig. 0Jd, and yield a n = 20.15 for 
the incident coherent state-a good agreement with the 
experimental value of 77, oxp = 4.02 |l9j . 
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FIG. 4: Reconstruction of the input-state photon-number dis- 
tribution: (a) Photon counting distribution from the NIST 
experiment ; (b) Reconstruction of the input state photon 
distribution (n « 20.15), assuming a quantum efficiency of 
0.2 and no dark counts. 



V. DARK COUNTS 

In addition to the finite quantum efficiency, the dark 
counts further make the results extracted from the exper- 
iment obscure. For a constant dark-count rate, however, 
a generalization of the inversion may be handled analyt- 
ically. Suppose in general that the probability of d dark 
counts is D(d), a discrete probability distribution which 
satisfies the normalization condition 



X>(d) = i 



(24) 



d=0 



We assume that the distribution D(d) represents proba- 
bilities of dark counts in a fixed interval of time r (du- 
ration of a single photon-counting experiment), and that 
these probabilities are independent of all other relevant 
observables such as incoming photon number, number 
of detected photons, etc. Therefore in the presence of 
dark counts, the conditional probability of detecting k 
photons, given n photons in the incident field Eq. @, is 
modified accordingly: 



P D (k\n) 



k 

£ 

d=0 



D(d) 



k-d 



v k - d (i-v) 



n—k-\-d 



(25) 



which expresses the probability Po(fc|n) as a sum (over 
all possible d) of the probability of registering k — d pho- 
tons from the incident field via true detection while the 
remaining d are registered via dark-count events. 

As long as the physical time scale for the emergence of 
a dark-count event is much smaller than r, and the events 
arrive at a constant rate A, the probability distribution 
D(d) (of the arrival of d dark-count events in the fixed 
time interval t) has a Poissonian-distribution limit: 



D(d) 



(26) 



where the average number of dark count is d = At. Since 
after interchanging the summation index d with k — d we 
can rewrite Eq. I|25|) as 



P D (k\n) = p^D(k-d)(^j V d (l 



7i — d 



(27) 



where the summation over d is in the form of a matrix 
multiplication [Eq. 



(-Po)fcn — £ ^kd Pdn, 



or in matrix notation, 

Pd =DP , 

where (Pn)^ = P D (k\n), and 

D kd ee D(k-d) 

-xr {^) k - d 
(k-d)\ 



(28) 



(29) 



(30) 



Throughout the paper, we adopt the convention that k\ = 
r(fc + 1) = oo for negative integer k, and correspondingly 
(^) = whenever b > a. In Eq. (|30|l . this convention 
implies that Dkd — for d > k. 

Now by inverting Eq. i[27J) in a matrix form, from the 
fundamental equation, 



P(k) = Y,PD(k\n)S(n) , 



(31) 



we can reconstruct the photon number distribution S(n) 
of the incident field given the vector of observed photon- 
count probabilities P{k), whose relationship to S(n) 
[given by Eq. I|31() ]. can again be written in matrix form 
as P = Pd-S, just as in Eqs. (20)-(21). The inverse of 
the coefficient matrix Po(k\n) = (Po)kn is found from 
Eq. 123), 



p-i D i 



(32) 



Let us first consider the inverse of the matrix D defined 
as in Eq. (30). It is not difficult to prove that the inverse 
of D is given by 



(i-j)l 



i-j 



(33) 



(i-i)l 

Combining Eq. (33) with Eq. (32) and Eq. (23), we find 



_ 1)k +n e Xr ^ Q) 



(i 



(l-77) fc (Ar)« f^ o (i-n)! 



i — n)\ 
V 



(34) 
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The infinite power series over i in Eq. (34) is absolutely 
convergent for all values of the argument 



(1 ~ V)^r 



(35) 



and can be expressed in terms of the generalized hyper- 
geometric function \F\ (a, 6, z) (which is an entire function 
of the argument z for all parameter values a, 6), giving 
a compact expression for the inverse matrix elements in 
the form: 



(Pn-'hn = (-1) 



k-\-n 



(l-T?) fc (Ar) 



n 



fc! (n-fc)! z " iFi(n + 1, n - k + 1, z) if fc sc n 
z fc iFi(fc + 1, k - n + 1, z) \ik>n 



W- 



(36) 



The inversion of Eq. (32) can now be effected by substi- 
tuting Eqs. (35)-(36) into the matrix formula 



(37) 



This suggests that the reconstruction of the photon 
statistics of the source may be possible with any given 
quantum efficiency of the detector, even in presence of 
dark counts [3(| • However, as the efficiency falls and the 
dark-count rate grows, more and more data points will 
be required for the inversion. 



VI. SUMMARY 

We have shown a way of quantifying the capability 
of photon-number discriminating detectors by the condi- 
tional probability, Q(n\k), given in Eq. ©. This condi- 
tional probability, however, can be obtained only when 
a priori knowledge of the input is available. Without 
a prior knowledge of the input state the results of pho- 
todetection can only eliminate the possibility of having 
certain numbers of photons in the source rather than de- 
termine them in the absence of dark counts. Hence, the 
term non-photon number- discriminating detector corre- 
sponds to the situation where the input state is unknown. 



On the other hand, the same photodetector could be a 
photon number- discriminating detector if the input state 
is known. Of course, the number resolving capability, 
without a prior knowledge of the input, can be measured 
qualitatively using the conditional probability P{k\n) of 
Eq. J5J) , which is given by the characteristics of the detec- 
tor itself. In particular, Ck = P(k\K) may serve the mea- 
sure of that capability. For the detector model considered 
in this paper, the number-resolving capability up to two 
photons can be written as Ck=2 = ?7 2 , and for Ci to be 
bigger than 0.5, r\ should be bigger than l/\/2 « 0.71. 
Similarly, Prj(k\k) of Eq. i|27|> plays the same role in the 
presence of dark counts. For example, with a given dark- 
count rate At = 0.5, the detector efficiency r\ should be 
more than 0.78 for Pn(k\k) to be larger than one half. 

We then introduced a simple matrix inversion formula 
in order to reconstruct the input-state photon distribu- 
tion from the measured photon counting distribution. 
Due to the form of the inverse matrix and its high sensi- 
tivity to the fluctuation, restrictions need to be imposed 
on the use of this inversion formalism. Alternatively, we 
can use a genetic algorithm that finds the distribution 
with maximum entropy after passing the y 2 goodness-of- 
fit test. We have shown good results for the reconstruc- 
tion of the input distribution using the data from a recent 
experiment at NIST. Finally, we have considered the re- 
construction task in presence of dark counts and derived 
a completely analytic formula for the inversion. A de- 
tailed analysis of the conditions for the inversion formula 
to be applicable will be discussed in future work. 
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